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THE e - REVISED SYSTEM OF THE RIGID BODY WITH 
THREE LINEAR CONTROLS 

Dan COMANESCU, Mihai IVAN and Gheorghe IVAN 

Abstract. In this paper we introduce the e - revised system associated to a Hamil- 
ton - Poisson system. The e - revised system of the rigid body with three linear controls 
is defined and some of its geometrical and dynamical properties are investigated. 

1 Introduction 

It is well known that many dynamical systems can be formulated using a Poisson 
structure (see for instance, R. Abraham and J. E. Marsden [1] and M. Puta [11]). 

The metriplectic systems was introduced by P. J. Morrison in the paper [8]. These 
systems combine both the conservative and dissipative systems. 

A metriplectic system is a differential system of the form x = PdH + gdC, where 
P is a Poisson tensor on a manifold M, g is a symmetric tensor of type (2, 0) on M, 
and H and C are two smooth functions on M with the additional requirements: 
(a) PdC = 0; (b) gdH = and (c) dC ■ gdC < 0. 

The differential systems of the form x = PdH + gdC which satisfies only the 
conditions (a) and (b) are called almost metriplectic systems ( see Fish, [2]; Marsden, 
[7] ; Ortega and Planas - Bielsa, [9] ) . An interesting class of almost metriplectic systems 
are so-called the revised dynamical systems associated to Hamilton- Poisson systems (see 
Gh. Ivan and D. Opri§, [5]). 

The control of the rotation rigid body is one of the problems with a large practical 
applicability. For this reason, in this paper we study the e - revised dynamical system 
associated to the rigid body with three linear controls. 

2 Almost metriplectic systems 

We start this section with the presentation of the concept of almost metriplectic 
manifold (see Ortega and Planas- Bielsa, [9]). 

Let M be a smooth manifold of dimension n and let C°°{M) be the ring of smooth 
real- valued functions on M. 

A Leibniz manifold is a pair (M, [•,•]), where [•, •] is a Leibniz bracket on M, that is 
[-, •] : C°°(M) x C°°(M) -> C°°(M) is a R - bilinear operation satisfying the following 
two conditions: 
(i) the left Leibniz rule: 

[fl-f2,h) = [fl,h)-h + fl-[f2,h] for all / 1 ,/ 2 ,/ 3 GC 00 (M); 
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(ii) the right Leibniz rule: 

[/i,/2-/ 3 ] = [/i,/2]-/3 + /2-[/i,/3] for all f 3 ,f 3 ,f 3 eC°°(M) 

where " •" denotes the ordinary multiplication of functions. 

Let P and g be two tensor fields of type (2, 0) on M and e G R be a parameter. 
We define the map [-, (-, -)] e : C°°(M) x (C°°(M) x C°°(M)) -> C°°(M) by: 

/ia)] e = P(df,dh 1 ) + eg(df,dh 2 ), for all f,h 1 ,h 2 eC QO (M). (1) 

Proposition 2.1. T/ie map [•,(•, -)]e given by (1) satisfy the following relations: 

(i) [afi + 6/ 2 , (hi, /i2)] e = a[fi, (hi,h 2 )} £ + &L/2, (hi, h 2 )] £ ; 

(ii) [f, a(h u h 2 ) + b(h[, h' 2 )] £ = a[f, fa, h 2 )] £ + b[f, (h[, h' 2 )} £ ; 
{in) [fh, (hi,h 2 )} £ = /[A, {h u h 2 )]e + (/ii, ^2)] £ ; 

(iv) [f, h(h u h 2 )} £ = h[f, {h u h 2 )]e + hiP(df, dh) + eh 2 g(df, dh), 
for all f, /1, f 2 , h u h 2 , h[, h' 2 G C°°(M) and a,kR. 

Proof. Applying the properties of the differential of functions and using that P 
and g are R- bilinear maps, it is easy to establish the relations (i) — (iv). □ 

We consider the map [[-, -]] e : C°°(M) x C°°(M) -» C°°(M) defined by: 

[[/, h}} £ = [f, (h, h)} £ , for all f,h€ C°°(M). (2) 

Therefore, the map [[•, -]] £ : C°°(M) x C°°(M) C°°(M) is given by: 

[[f,h]] £ = P(df,dh) + eg(df,dh), for all f,heC°°(M). (3) 

Proposition 2.2. T/ie bracket [[■, -]] e on M ozwen 6j/ (3) verify the right Leibniz rule: 
[[f,hh% = h[[f,h% + h'[[f,h]] £ , for all f,h,ti GC°°(M). 

Proof. Indeed, [[/, = [f,(hh',hh% = [f,h(h',h%. Putting h x = h 2 = h! in 

the relation (in) from Proposition 2.1, we have succesive: 

[/, h(h>, h% = h[f, (h f , h')} + h'P(df, dh) +eh'g(df, dh) = h[f, (h', h% + h'(P(df, dh) + 
eg(df, dh)) = h[f, (h', h% + h'[f, (h, h)} £ = h[[f, h% + h'[[f, h}} £ . □ 

By Proposition 2.1. (i), (ii) and (hi) and Proposition 2.2, we have that the map 
[[•, -]] e given by (3) is a Leibniz bracket on M. Hence, [[•, -]] e defines a Leibniz structure 
on the manifold M and (M, P, g, [[•, -]] e ) is a Leibniz manifold for each e G R. 
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A Leibniz manifold (M, P, g, [[■, -]] £ ) such that P is a skewsymmetric tensor field 
and g is a symmetric tensor field is called almost metriplectic manifold. In other words, 
given a skewsymmetric tensor field P of type (2, 0) and a symmetric tensor field g of 
type (2, 0) on a manifold M, we can define an almost metriplectic structure on M. 

If the tensor field P is Poisson and the tensor field g is nondegenerate, then 
(M, P, g, [[-,-]] e ) is a metriplectic manifold, see Ortega k, Planas - Bielsa [9]. 

Proposition 2.3. Let (M,P,g, [[■, -]] £ ) be an almost metriplectic manifold. If there 
exist h 1 ,h 2 e C°°(M) such that P(df, dh 2 ) = and g(df, dh x ) = for all f G C°°(M), 
then the bracket [[■, •]] given by (3) satisfies the relation: 

[[/, h! + h 2 }] £ = [f, (hi,h 2 )\ £ , for all f G C°°(M). (4) 

Proof. Indeed, [[f,h 1 + h 2 ]] £ = P(df,d(h 1 + h 2 )) + £g(df,d(h 1 +h 2 )) = P(df,dh 1 + 
dh 2 )+eg(df,dhi+dh 2 ) = P(df,dhi)+P(df,dh 2 )+eg(df,dhi)+eg(df,dh 2 ) = P(df,dhi)+ 
eg(df,dh 2 ) = [f,(h l ,h 2 )} £ . □ 

Let (M, P, g, [[•, -]] £ ) be an almost metriplectic manifold and let hi,h 2 € C°°(M) 
two functions such that P(df,dh 2 ) = and g(df,dhi) = for all / G C°°(M). The 
vector field Xh 1 h 2 given by: 

X hlh2 (f) = [[f,h 1+ h 2 ]] £ for any / G C°°(M) 

is called the Leibniz vector field associated to the triple (hi, h 2 ,e) on M. 
Taking account into Proposition 2.3 and (1), Xh 1 h 2 1S given by: 

X hlh2 (f) = [f,(hi,h 2 )} £ = P(df,dh 1 ) + eg(df,dh 2 ), for all / G C°°(M). (5) 

In local coordinates on M, the differential system given by: 

x l = [[x l ,h l +h 2 ]] £ = [x i ,(h l ,h 2 )] £ (6) 

where 

[x\(h l ,h 2 )] £ = X hlh2 (x*) = P^+eg^, i,j = l^ (7) 

with P % i = P(dx l , dxi) and g i: > = g(dx l , dx^), is called the almost metriplectic system 
on M associated to the Leibniz vector field X^^ with the bracket [[•, -}] £ . 

We denote the matrix of the tensor fields P and g respectively by P = (P l i) and 
9 = (q 13 )- We have that P is a skewsymmetric matrix and g is a symmetric matrix. 

We give now a way for to produce almost metriplectic manifolds. 

Proposition 2.4. For a skewsymmetric tensor P of type (2, 0) on a manifold M 
and two functions h x ,h 2 G C°°(M) such that P(df,dh 2 ) = for all f G C°°(M), 
there exists a symmetric tensor g of type (2,0) on M such that g(df,dh\) = for all 
f G C°°(M) and (M,P,g, [[■, -]] e ) is an almost metriplectic manifold. 
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Proof. In a system of local coordinates on M, let g = (g lJ ) the matrix of the 

symmetric tensor g which must to be determined. Then, the components g lJ , i,j = l.n 

Oh 

verify the system of differential equations g lJ — = 0, i,j = l,n. 

dh 

In a chart U such that — — (x) / , the components g %3 are given by: 

oxi 

9 U (*) = ~ E Qf 

k=l, k^i ox (Q\ 
H( ^ dhidhi . ... y ) 

- 9 {X) = dx^dxl' fOT 1 + 3 
Applying now Proposition 2.3 we obtain the result. □ 

Proposition 2.4 is useful when we consider the e- revised system of a Hamilton- 
Poisson system. 

For this, let be a Hamilton-Poisson system on M described by the Poisson tensor 
P having the matrix P = (P iJ ) and by the Hamiltonian function h\ £ C°°(M) with 
the Casimir function h 2 G C°°(M) ( i.e. -P^f^j = for i,j = 1, n ). The differential 
equations of the Hamilton-Poisson system are the following: 

f)h 

<,J=I^- (9) 

Using (8), we determine the matrix g = (g tJ ) and we have: 

f)h 

.'/"— iJ = ^- (10) 

Applying now Proposition 2.4, for each e £ R, we obtain an almost metriplectic 
structure on M associated to system (9). The differential system associated to this 
structure is called the e - revised system of the Hamilton - Poisson system. 

Hence, the s - revised system of the Hamilton - Poisson system defined by (9) is: 

•i , ii^2 . . i / 1 ~t \ 

x * = p*3 +£g V l,] = l,n. (11) 

The terms fl^f^j, i,j = l,n from the e - revised system (11) describe a cube pertur- 
bation of the Hamilton - Poisson system. 



Remark 2.1. We observe that the 0- revised system (11) coincide with the Hamil- 
ton - Poisson system (9). □ 
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3 The £ - revised system associated to the rigid body with 
three linear controls 

The rigid body equations with three linear controls ( see, M. Puta and D. Comanescu 
[12] ) are given by: 

x 1 = (a 3 — a 2 )x 2 x 3 + cx 2 — bx 3 

x 2 = (ai — as)x 1 x 3 — cx 1 + ax 3 (12) 
x 3 = (ci2 — ai)x 1 x 2 + bx 1 — ax 2 

where x(t) = (x 1 (t) , x 2 (t) , x 3 (t)) G R 3 and a\ = -!-, a 2 = -!-, a 3 = -!- with I± > I 2 > 

h h h 

Is > ( I\,l2,l3 being the principal moments of inertia of the body ) and a, b, c € R 
are feedback parameters. We have < a± < a 2 < a^. 

The dynamics (12) is described by the Poisson tensor IT and by the Hamiltonian H 
on R 3 given by: 

/ -x 3 x 2 \ 
n(x) = i 3 -x 1 , (13) 
\ -x 2 x l / 

H(x) = ^[a^x 1 ) 2 + a 2 {x 2 ) 2 + a 3 (x 3 ) 2 ] + ax 1 + bx 2 + cx 3 . (14) 

Using (13) and (14), the dynamics (12) can be written in the matrix form: 

x(t) =n(x(t)) • VH(x(t)), (15) 

where x(t) = (x 1 (t) , x 2 (t) , x 3 (t)) T and VH(x(t)) is the gradient of the Hamiltonian 
function H with respect to the canonical metric on R 3 . 

Therefore, the dynamics (12) has the Hamilton-Poisson formulation (R 3 ,IT, H), 
where II and H are given by (13) and (14). 

The function C G C°°(R 3 ) given by: 

C(x) = \[(x 1 ) 2 + {x 2 ) 2 + (x 3 ) 2 ] (16) 

is a Casimir of the configuration (R 3 ,II), i.e. 

C(x) ■ VH{x) = O. (17) 

Applying the relations (8) for P = H,hi(x) = H(x) and h 2 (x) = C(x), the sym- 
metric tensor g is given by the matrix: 

-(a 2 x 2 + b) 2 — (a^x 3 + c) 2 (aix 1 + a)(a 2 x 2 + b) (aix 1 + a)(asx 3 + c) 

(aix 1 + a)(a 2 x 2 + b) -(aix 1 + a) 2 - (a 3 x 3 + c) 2 (a 2 x 2 + 6)(a 3 x 3 + c) 
(aix 1 + a)(a 3 x 3 + c) (a 2 x 2 + 6)(a 3 x 3 + c) -(aix 1 + a) 2 - (a 2 ^ 2 + b) 2 
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dhi i dhi 2 , dh\ ■> 

since — r = a\x + a, — = a 2 x + 6, — ? = a 3 x° + c. 
ax 1 ctar ax- 1 

We have 

a(x) • Vn 2 (x) = (v 1 (x),v 2 {x),v 3 {x)) T (18) 

where 

t>i(x) = -[(a 2 x 2 + 6) 2 + (a 3 x 3 + c) 2 ]x 1 + (aix 1 + a)[(a 2 x 2 + 6)x 2 + (a 3 x 3 + c)x 3 ] 
^(^c) = -[(aia; 1 + a) 2 + (a 3 x 3 + c) 2 ]x 2 + (a 2 x 2 + 6)[(aix 1 + a)x 1 + (a 3 x 3 + c)x 3 ] 
i^O^) = -[(ai^ 1 + a ) 2 + (a 2 x 2 + 6) 2 ]x 3 + (a 3 x 3 + ^[(aix 1 + a)x 1 + (a 2 x 2 + 6)x 2 ] 

(19) 

The e - revised system associated to dynamics (12) is: 

f x 1 = [(a 3 — a 2 )x 2 x 3 + cx 2 — 6x 3 ] + ev\{x) 

[(ai — a 3 )x 1 x 3 - cx 1 + ax 3 ] + ev 2 (x) (20) 
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x 

x 3 = [(a 2 — ai)x 1 x 2 + ftx 1 — ax 2 ] + ev 3 {x) 

The differential system (20) is called the e - revised system of the rigid body with 
three linear controls. Taking a = b = c = and e = 1 in (20), we obtain the revised 
system of the free rigid body, see [5]. 

Vector writing of the e- revised system (20). We introduce the following no- 
tations: 

x = (x\x 2 ,x 3 ), v = (vi,v 2 ,v 3 ), a=(a, b, c), m(x) = (aix 1 + a, a 2 x 2 + 6, a 3 x 3 + c). 

For all u = (u±, -u 2 , u 3 ), w = (i«i, n> 2 , W3) G R 3 > the following relation holds: 

u • [w x (u x w)] = (u x w) 2 , with (u x w) 2 = (u x w) • (u x w) (21) 
where " x " and " • " denote the cross product resp. inner product in R 3 ; that is: 

U X W = (li 2 U> 3 — U3W2, U3W1 — U1W3, U\W2 — U 2 W\), U • W = U\W\ + U2W2 + U3W3. 

With the above notations, the dynamics (12) has the vector form: 

x = x x m(x). (22) 

It is not hard to verify the following equality: 

v = (x x m(x)) x m(x). (23) 

Using the relations (22), (23) and (20), we can written the e- revised system (20) in 
the vector form: 

x = x x m(x) + e[(x x m(x)) x m(x)]. (24) 
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4 The equilibrium points of the e - revised system 

The equilibrium points of the Hamilton - Poisson system (12) ( or (22) ) are solu- 
tions of the vector equation: 

x x m(x) = 0. (25) 

The equilibrium points of the e - revised system (20) ( or (24) ) are solutions of the 
vector equation: 

x x m(x) + e[(x x m(x)) x m(x)] = 0. (26) 



Theorem 4.1. The Hamilton - Poisson system (12) and its revised system (20) 
have the same equilibrium points. 

Proof. Let xo be an equilibrium point of the system (12). According with (25) 
follows xo x m(xo) = 0. We have that xo is a solution of the vector equation (26), since 
xo x m(xo) + (xo x m(xo)) x m(xo) = + x m(xo) = 0. Hence xo is an equilibrium 
point of the e- revised system (20). 

Conversely, let xo be an equilibrium point for (20). Using (25) it follows 

(a) x x m(x ) + e[x x m(x )) x m(x )] = 
The relation (a) can be written in the form: 

(b) x x m(x ) - e[m(x ) x (x ) x m(x ))] = 
Multiplying the relation (b) with the vector xo, we obtain: 

(c) x • (x x m(x )) - ex • [m(x ) x (x ) x m(x ))] = 0. 
Using the equality (21), the relation (c) is equivalent with: 

(d) - e(x x m(x )) 2 = 0. 

From (d) (if e ^ 0), follows xo x m(xo) = 0, that is xo is an equilibrium point for (20). 
□ 

The equilibrium points of the dynamics (12) are well-known (see M. Puta and D. 
Comanescu, [12]) and these are presented in the following proposition. 

Proposition 4.1. ([12]) The equilibrium points of the Hamilton - Poisson system 
(12) are the following: 

(i) ei = (0,0,0); 

(ii) e 2 = (y~ — >\ 'T~~ — ) f or AeR\{oi,o 2) a 3 }; 

A — a\ A — a 2 A — a 3 

b c 

(Hi) e 3 = (a, , ) for a G R, if a = 0; 

a 2 — a\ a\ — a 3 

a c 

(iv) e4 = ( , a, ) for a G R, if 6 = 0; 

a 2 — ai a 3 — a 2 

(v) e 5 = ( — , ,a) for a G R, if c = 0. 

cl\ — a 2 a 3 - a 2 
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By Theorem 4.1, the equilibrium points of the e- revised system (20) are ei, e$ 
indicated in the Proposition 4.1. 

It is well-known that the dynamics (12) have the first integrals H and C given by 
(14) and (16). These first integrals may be written thus: 

H(x\x 2 ,x 3 ) = -x-I _1 x + a-x and C(x 1 , x 2 , x 3 ) = -x 2 (27) 

where I is inertia tensor and I -1 is its inverse. We have: 

dH dC . . , . 

— — (x) = m(x) • x and — (x) = x • x. (28) 

\Jbh HI/ 

dH 

Indeed, = (aix 1 + a)x 1 + fax 2 + b)x 2 + (a^x 3 + c)x 3 = m(x) • x and 

dC 1-1 2-2 s • s 

= X X + X X + X X = X • X. 

dt 

Theorem 4.2. (i) For each e € R, i/ie function H given by (14) is a /irst integral 
for the e- revised system (20). 

(ii) If x : R — > R 3 zs a solution of the e- revised system, then: 

|(ix 2 ) = -e(xxm(x)) 2 . (29) 

(Hi) For e G R*, the function C is not a first integral for the e- revised system. 

Proof, (i) Multiplying the relation (24) with the vector m(x), we have: 

m(x) • x = m(x) • (x x m(x)) + em(x) • [(x x m(x)) x m(x)] = 0. Applying now (28), 
dH 

we obtain — — = m(x) • x = 0. Hence H is a first integral for the system (20). 

(ii) Multiplying the relation (24) with the vector x, we have 
x • x = x • (x x m(x)) + ex • [(x x m(x)) x m(x)] = —ex • [m(x) x (x x m(x))]. 

Using now the equality (21), we obtain x • x = — e(x x m(x)) 2 . Then, we have 

4-(-x 2 ) = x • x = -(x x m(x)) 2 . 

dV2 1 \ \ >) 

(Hi) This assertion follows from the second relation of (28) and (ii). □ 
Remark 4.1. The function H given by (14) can be put in the equivalent form: 

H(x 1 ,x 2 ,x 3 ) = l[a 1 (x 1 + -) 2 + a 2 (x 2 + -) 2 + a 3 (x 3 + -) 2 ]-l(- + - + -). (30) 
z a\ 02 03 z ai 02 (13 

For a given constant k G R, the geometrical image of the surface: 
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H(x\x 2 ,x 3 ) = k 

is an ellipsoid, since a± > 0, 02 > 0, 03 > 0. □ 

Proposition 4.2. The set of equilibrium points which belong to the ellipsoid 
H(x 1 , x 2 , x 3 ) = k, is finite. 

Proof. Following the description of the equilibrium points given in Proposition 4.1, 
we remark that: 

(i) the equilibrium points of the form e^ ( similarly, for e^ and e^ ) make a 
straight line; the intersection between a straight line and an ellipsoid have at most two 
points; we deduce that on the chosen ellipsoid there exist at most two points of the 
form e^. 

(ii) the equilibrium points of the form e2 can be obtained by solving with respect 
A the following equation: 

2 A — a\ A — a>2 A — 03 A — a± A — 02 A — 03 

The above equation is equivalent with the determination of roots of a polynomial 
of degree at most 6; therefore on the chosen ellipsoid there exist at most 6 equilibrium 
points of the form e<i . □ 



5 The behaviour of the solutions of the e - revised system 

Theorem 5.1 (i) The solutions of the e-revised system are bounded. 

(ii) The maximal solutions of the e-revised system are globally solutions (i.e. these are 

defined on R). 

Proof (i) Given a solution of (20), there exists a constant k such that its trajectory 
lie on the ellipsoid H(x 1 , x 2 , x 3 ) = k. From this we deduce that all solutions are 
bounded. 

(ii) Let x : (m, M) C R — > R 3 be a maximal solution. We assume that x is not 
globally. It follows m > — 00 or M< 00. In these situations, we known that there 
exists k £ R such that H(x 1 , x 2 , x 3 ) = k for all t 6 R and the graph of the solution is 
contained in a compact domain. According with [6] (theorem 3.2.5, p. 141) we obtain a 
contradiction with the fact that x admit a prolongation on the right or the left (also, 
can be applied the theorem of Chilingworth (1976), see theorem 1.0.3, p. 7 in [3]). □ 

In the sequel we study the asymptotic behaviour of the globally solutions of the 
e-revised system. 
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Denote by E the set of equilibrium points of the e-revised system (20) and by T the 
trajectory of a solution x : R — > R 3 of (20). By theory of differential equations (see 
[10] p. 174-176), the w-limit set and a-limit set of T are: 

u(F) = {y £ R 3 / 3t n — > oo such that x(i n ) -> y}, 

a(r) = {z E R 3 / 3i n -> -co such that x(i„) -> z}. 

Theorem 5.2 Let x : R — > R 3 be a solution of the e-revised system with e ^ 
0. There exist the equilibrium points x m ,XM G R 3 of the system (20) such that 
lim^-oo x(i) = x M and lim^^ x(i) = x m . 

Proof The theorem is proved in the following steps: 

(i) a (r) ^ and co(F) + 0. 

(ii) a(r)nw(r)cE. 

(hi) The sets a(F) and u>(F) contains exactly one element. 

Taking account into that each solution is bounded (hence it is contained in a com- 
pact domain) and applying theorem 1, p. 175 in [10], we obtain immediately the 
assertions (i). 

(ii) For demonstration consider the case when e > 0. Using the relation (29), we 
deduce that the function t — ► x 2 (i) is a strictly decreasing function. Being bounded it 
follows that there exists lim^ 00 x 2 (t) = L and L is hnite. 

For each y G to(F) there exists the sequence t n — > oo such that x(i„) — ► y. Then 
x2 (i«) — ^ y 2 and hence y 2 = L. 

By theorem 2, p. 176 in [10], we have that the trajectory Ty of the solution xy 
which verifies the initial condition xy(0) = y, satisfies the relation Ty C u>(T). 

If we assume that y is not an equilibrium point, then we deduce (using the relation 
(29)) that for t > we have Xy(t) < L and this is in contradiction with the above 
result. Therefore, we have uj(T) C E. 

Similarly, we prove that a(T) C E. Hence the assertion (ii) holds. 

The case e < is similar. 

(iii) There exists a constant k such that the sets a(T) and o>(r) are included in the 
ellipsoid H(x 1 ,x 2 ,x 3 ) = k. By (ii), we deduce that a(F) and u)(F) are included in the 
set of equilibrium points which lies of the above ellipsoid. On the other hand, applying 
Proposition 4.2 and using the fact that a(r) and u>(T) are connected (see theorem 1, 
p. 175 in [10]), we obtain that a(T) and uj(F) are formed by only one element. □ 

Remark 5.1 Using the relation (29) it is easy to observe that the following asser- 
tions hold: 

(i) if £ > then x 2 M > x^; 

(ii) if e < then x 2 M < x^. □ 
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As an immediate consequence we obtain the following theorem. 

Theorem 5.3 If e ^ 0, then for each solution x : R — > R 3 of the e-revised system 
we have: 

if t^oo => d(x(t),E)-»0 , . 

if t-» -oo =► d(x(t),E) -v0. 1 j 

□ 

Remark 5.2 From Theorem 5.3 follows that the set E of the equilibrium points 
is an attracting set (see definition 2, p. 178 in [10]) and also is a reppeling set (see [3], 
p. 34). Thus, the space R 3 is simultaneously a domain of attraction and a domain of 
repulsion of E. □ 

6 The Lyapunov stability of equilibrium points of the e - 
revised system in the case e > 

The stability of the point e\ = (0,0,0). We have the following results. 
Theorem 6.1 The equilibrium point e\ is Lyapunov stable. 

Proof Let 7 > 0, to G R and xo G R 3 such that |xo| < 7, where | • | denotes the 
euclidian norm in R 3 . Denote by t — > x(t,to> x o) the solution of the e-revised system 
which verifies the initial condition x(0, to> x o) = x o- 

Using the relation |x(i,io, x o) = A/x 2 (i, to, xo) and according with the relation 
(29), we observe that the function t — > x(t, io) x o) is a decreasing function and hence 
we have: 

|x(t, t , x o)| < l x o| < 7 for t > to- 
Then (see [4], p. 22) we have that e\ is a Lyapunov stable equilibrium point. □ 

Remark 6.1 The equilibrium point e\ is not asymptotical stable. 
Indeed, ifa = 6 = c = then the coordinates axis are formed from equilibrium 
points. If al least of one of the numbers a, b, c is non null, then: 

if |A|-oc => (t-^— >T^— >T^— ) ^ (0 '°' 0) - 
A — a\ A — ai A — C13 

Hence, in all neighbourhood of e\ there exist an infinity of equilibrium points. □ 

The stability of the point xo = (— — , — — , — — ). The equilibrium point xo is 

a\ (12 0.3 

an equilibrium point of the form e2 and it is obtained for A = 0. 
Theorem 6.2 The equilibrium point xq is Lyapunov stable. 
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Proof Using the relation (30) and the inequality < a\ < a-i < as, we deduce: 

y|x-X^| < H(X.) -Hfa) < y|x-5c^| 

For to G R and xo G R 3 denote with x(i, *o,xo) a solution of e-revised system which 
verifies the initial condition x(0, io> x o) = x o- 

Let 7 > and <5( 7 ) = f^. Let xo G R 3 such that: 

ff(xo) - fl"(x5) < <5( 7 ). 

From the fact that i? is a first integral we deduce that: 

fT(x(i, t , x ) - fT(xo) = # (x ) - ff(xo). 

Hence for all t G R the following inequality holds: 

y|x-xo"| < 5(7) 

and we obtain that xo is Lyapunov stable. □ 

Remark 6.2 The stable equilibrium point xo realizes the absolute minimum of the 
function H. □ 

The unstability of equilibrium points of the form e2 with A € (0, a±). For 
the demonstration of this results we use the Theorem 6.3 and Lemma 6.1. 

Theorem 6.3 If xo £ E such that there exists y G E with the properties: 
(i) H(y) = H(x ) and (ii) |y| < |x 
then xo is an unstable equilibrium point. 

Proof For k G R denote by E& = {x G E / H(x) = k}. The set E^( Xo ) is finite (by 
Proposition 4.2). We denote: 

70 = min{|x - x |/x G E h{Xq) - {x }} 

Let z G R 3 such that H(z) = i^(xo) and |z| < |xo|. Then: 

lim x(t, 0,z) G E 

t—>oo 

and 

if t > => |x(t,0,z)| < |z| 
and we deduce that there exists t z > such that: 

|x(t,0,z) -x | > y if t > t g 
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It follows that xo is unstable. □ 

We assume that (a, b, c) ^ (0, 0, 0) and we introduce the notation: 

, a b c . , , 

e 2 \ = (t , t , t ) for all A € R — {ai,a 2 ,a 3 } 

A — a± A — Ci 2 A — (2 3 



Lemma 6.1 (i) If a < \i < a±, then \e 2a \ < l^l- 
(ii) If a,n > 0, (^) 2 > ^ and #(e 2(T ) = H(e 2 ^), then |e 2(J | > |e 2M |. 
(hi) If < cr < ai < a 3 < /i and H(e 2a ) = H(e 2fl ), then |e 2cr | > |e 2jU |. 

Proof (i) Consider the function g : (— oo, ai) — > R given by: 

The derivative of the function g is: 

2a 2 2b 2 2c 2 



9'W = - 



(A-ai)3 (A-a 2 )3 (A - a 3 ) J 



We observe that g'(X) > and we obtain that g is a strictly increasing function. We 
have: 

g(<r) = \e 2(7 \ 2 , g(ji) = \e 2fl \ 2 

and we obtain the desired result. 

(ii) From hypothesis H(e 2a ) = H(e 2fl ) follows that there exists a constant q > 
with the following properties: 

la 2 lb 2 1 c 2 q 

+ —7 + 



oi (cr - ai) 2 a 2 (cr - a 2 ) 2 a 3 (cr - a 3 ) 2 cr 2 
la 2 lb 2 1 c 2 a 

+ —7 ^ + 



ai{jx-ai) 2 a 2 {jx-a 2 ) 2 a 3 (/i - a 3 ) 2 ^ 2 
Using a\ < a 2 < a 3 , we obtain the inequalities: 

I |2 ^ °19 | |2 ^ «3<? 

e2<r > , e 2At < 

o n 

and we observe that the assertion (ii) holds. 

(hi) This assertion follows immediately from (ii). □ 

Theorem 6.4 The equilibrium point e 2 \ with < A < a\ is unstable. 
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Proof Consider the function h : (—00, a±) |J( a 3> 00 ) ~~ y R- given by: 



^(a) = #(e 2(T ). 



Using the relation (30) for H, we find: 



a 



2 




2 ai(<7 — ai) 2 



+ 



a 2 (a - a 2 ) 2 



+ 



a 3 (cr-a 3 ) 2 2 ai a 2 03 



The function h have the following properties: 

• is an absolute minimum point. 

• lim^-oo h(a) = lim^oo h(a) = 0. 

• lim CT ^ ai h(a) = lim CT _» a3 h(a) = 00. 

• h is strictly decreasing on (— 00, 0), strictly increasing on (0, a\) and strictly de- 
creasing on (03,00). 

The demonstrations divided on three cases. 

(I) Assume that h(X) < 0. In this situation there exists a < < A < a\ such that 
h(X) = h(a) and imply H{e,2\) = i7(e 2CT ). Hence the equilibrium points e 2 >, and e 2cr 
belong to same ellipsoid. 

On the other hand, by Lemma 6.1 (ii), follows |e 2(T | < |e 2 ^|. Applying now Theorem 
6.3, deduce that e 2 A is an unstable equilibrium point. 

(II) Assume that h(\) = we have H{e2\) = H(0, 0, 0) and it is clearly that 
|(0,0,0)| < |e2A I - By Theorem 6.3 we find the desired result. 

(III) Assume that h(X) > 0. Then there exists a > a 3 such that h(\) = h(a) and 
hence H(e.2\) = iJ(e 2(T ). Applying Lemma 6.1 (iii) follows |e 2 ^| > |e 2(J | and by Theorem 
6.3 we deduce that e 2 A is unstable. □ 

The stability of equilibrium points of the form e 2 with A < 0. 

Theorem 6.5 The equilibrium points of the form e 2 with A < are Lyapunov 
stables. 

Proof Let A < and the equilibrium point xo = (jz^, jz^t X=oj) °^ f° rm e 2- 
It is well-known that the study of stability of xo in the Lyapunov sense is equivalent with 
the study of stability of the null solution (0, 0, 0) for the differential system obtained 
from the e-revised system by transformation of variables: 



z = 



x - x 



The system obtained in this manner is called the perturbed e-revised system. 
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Consider the function K : R 3 — > R given by 

K(z) = -z-I- 1 z--z 2 

Since the tensor I" 1 is strictly positive definite and A < we obtain that K is a 
quadratic form strictly positive definite. 

Next we prove that if z : R — > R 3 is a solution for the perturbed ^-revised system, 
then: 

><» 

By a direct computation and taking account into the relations (27) we have: 

K(z) = H(x) - AC(x) - ^x • I _1 x - a • x + -x 2 , 

Applying now Theorem 4.2, we obtain: 

|K(z(t))=eA(xxm(x)) 2 

and follows that ^K(z(t)) < 0, since e > 0, A < 0. 

It is easy to see that K*(t) = K(z(t)) is a strictly decreasing function. By theorem 
1.1, p. 21 in the paper [4] we deduce that xo is Lyapunov stable. □ 

Conclusion - the stability problem for the Hamilton Poisson system (12) 
versus the e- revised system (20) with e > 

Concerning to the equilibrium points of the system (12) are established the following 
results (see, theorem 1.1, [12]): 

(1) ei is Lyapunov stable; 

(2) e2 are Lyapunov stables for A G (— oo, a±) U (03, 00); 

(3) e3 are Lyapunov stables; 

(4) e4 are unstables; 

(5) es are Lyapunov stables. 

By Remark 4.1 (see, [12]), there exist cases for which the problem to decide the 
nonlinear stability or unstability are not discussed. 

For the stability of equilibrium points of the e- revised system (20) with e > have 
proved the following assertions: 

(1) ei is Lyapunov stable; 

(2) the equilibrium points of the form e2 with A < are Lyapunov stables; 

(3) the equilibrium points of the form e2 with < A < a\ are unstables. 
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